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ABSTRACT 

The Large Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) is a dedi- 
cated spectroscopic survey telescope being built in China, with an effective aperture 
of 4 meters and equiped with 4000 fibers. Using the LAMOST telescope, one could 
make redshift survey of the large scale structure (LSS). The baryon acoustic oscillation 
(BAO) features in the LSS power spectrum provide standard rulers for measuring dark 
energy and other cosmological parameters. In this paper we investigate the meaure- 
ment precision achievable for a few possible surveys: (1) a magnitude limited survey 
of all galaxies, (2) a survey of color selected red luminous galaxies (LRG), and (3) a 
magnitude limited, high density survey of 2; < 2 quasars. For each survey, we use the 
halo model to estimate the bias of the sample, and calculate the effective volume. We 
then use the Fisher matrix method to forecast the error on the dark energy equation 
of state and other cosmological parameters for different survey parameters. In a few 
cases we also use the Markov Chain Monte Carlo (MCMC) method to make the same 
forecast as a comparison. The fiber time required for each of these surveys is also 
estimated. These results would be useful in designing the surveys for LAMOST. 

Key words: large scale structure; cosmological parameters. 



1 INTRODUCTION 



Baryon acoustic oscillations (BAO) before recombination 
left wiggling features in the matter power spectrum. At 
large scale where the evolution of the density perturba- 
tion is still linea r, such features is pre s erved in the galaxy 
power spectrum ( Eisenstein fc Hu|[l998l : iMeiksin et. all 19991 : 
lEisenstein et afl 



l2005h an d the Sloan D igital 
vey datallEisensteinI I200F 



Sk y Surv ey 
Hut"siTl2005al : 



(SPS S) 



Huetsi I l2005bl: 



19971 ) 



This provides us with a well cal- 
ibrated standard ruler, which enables precise measure- 



Darameters (lEisenstein & Whitd 2004: EisensteinI 2005a: 


Sec & EisensteinI 20031, 20071: 


Koehler et all 2007: iLinder 1 


20031: Guzik & BernsteinI 2007 


: Blake & Glazebrookl 20031: 


Blake et al.l I2OO61). Like cosmic microwave backEround 



(CMB), the physics involved is relatively clean and well un- 
derstood, hence it is arguably less affected by unknown sys- 
tematic errors, which might undermine empirical-rule based 
methods such as the type la supernovae (SNIa) and the clus- 
ter abundance measurements. 

Recently, such a BAO feature have bee n observed 
in the 2dF galaxy redshift survey (2dFGRS) (|Cole et all 



'Tegmark ct al.' '2006': 'Percival et al. H2007l : lOkumura et all 
^2008; C abre fcGaztanaga 2008a J). A number of more 
powerful BAO s urveys are being plann ed. Some exam- 
ples are WiggleZ (|Glazebrook et aI.!l2 007^ SDSS -3 (BOSS) 
n, HETDEX 0, WFMOS (B assctt eTal] |2005| ). ADEPT 
0, Space/Euclid (Cima tti et all I200i') in the optical/IR 
and FAStB HSHStPeterson. Bandura fc PenllToOQ ). SKA 
l|Abdalla fc RawUngjiioOSl ) Fl in the radio. 

In this paper, we make forecasts on cosmological con- 
straints from BAO measurement in future surveys with 
the Large Sky Area Mul ti-Object Sp ectroscopic Telescope, or 
LAMOSTFI ((Chuiil998i : iStondl2008l ). a telescope being built 



http: //cosmology. lbl.gov/BOSS/\ \http://u 
http://www.as.utexas.edu/hetdex/ 



Ksdss3.org/\ 



http: //www7.nationalacademies. org /ssb/be^nov^ 2006.bennett.pdf\ 
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in China. The LAMOST is a 4 meter Schmidt telescope 
with a field of view of 20 deg^ Equipped with 4000 opti- 
cal fibers that are individually positioned by computer dur- 
ing observation, it is able to take spectra of 4000 targets 
simultanously. More details of the LAMOST telescope are 
described in Appendix A. The LAMOST is ideally suited 
for conducting large scale redshift surveys. 

We expect that once LAMOST becomes operational, it 
will carry out several different survey projects. The design 
of such surveys should maximize the scientific output for a 
given amount of observing fiber time. Some obvious choices 
include a magnitude limited general galaxy redshift survey, 
and a magnitude limited quasar survey. For the purpose of 
BAO measurements, it is desirable to have additional sur- 
veys for targets at relatively higher redshift. Often found to 
be the bright central galaxies of galaxy groups and clusters, 
the luminous red galaxies (LRG) ([E iscristcin ot al. 2001j; 
iBrown et alll2003l : IZehavi et al.ll2005ah can be detected at 
higher redshifts than a typical galaxy for a given apparent 
magnitude. Quasars can be observed at even higher red- 
shifts. In the SDSS surveys, quasars were v ery sparse, mak- 
ing it diffic ult to use in BAO measurement ([Schneider et al.l 
I2OO5I . I2OO7I ). However, it is expected that the density of ob- 
servable quasars with the LAMOST is much higher, since 
a lower limiting flux can be reached. Therefore, we also 
consider quasar samples in our study, with quasar them- 
selves as the tracer of the large-scale structure (LSS). The 
absorption lines (Lya fo rest) in the quasar spectra also pro- 
vide useful information (IWhite ll2003l : lMcDonald et al.ll2006l : 
iMcDonald fc Eisensteinll2007l ). which will be considered in 
our future work. We describe these surveys in more detail 
(including the selection of the sample, the required observ- 
ing fiber hours, and the estimate of bias of the sample) in 
Appendix B. 

We will make forecasts of the measurement errors on 
the dark energy equation of state and other cosmological pa- 
rameters with these surveys, and also estimate the resources 
required for these surveys. We w i ll first use the Fisher matrix 
method l|Seo fc Eisensteinll2003l . l2007l ). Then for a few cases 
we also use a Markov Chain Monte Carlo (MCMC) sim- 
ulation to make the forecast. The advantages of the Fisher 
matrix method are that it can easily be used to explore large 
parameter space and that the relationship between the in- 
put and output is very clear. The MCMC method does not 
assume the likelihood to be Gaussian and thus can probe 
the full shape of the likelihood surface. 

Previous studies have shown that it is feasible with 
the LAMOST to obtain better precisions in cosmolog- 
ical parameters than ongoing surveys like t h e SDSS 
jFeng, Chu. fc Yand I2OO0I: ISun. Su fc Fanl I2OO6I : iLi et all 
l200a ) The studv bv iFeng. Chu. fc Yand (120001 ) were con- 
ducted before the importance of B AO in d ark energy mea- 
surement was widely recognized. ISun. Su . & Fan (2006) 
consi dered the Alcock -Paczynski test (|Alcock fc Paczvnskii 
I1979I ) in real space (|Matsubara fc Szalavl I2OO3I ) . with a 
method different from the one we use. In all of these studies, 
a single survey of galaxies up to _B = 20.5 were considered, 
assuming an total survey area of 15000 deg'^, a total number 
of galaxies 10^, and a galaxy sample bias of 1. 

In this paper, we try to make a realistic assessment of 
LAMOST surveys. We consider different samples that could 
be collected by the LAMOST telescope within a reasonable 



amount of observation time. The surveyed sky area is as- 
sumed to be about 8000 deg^, for which the SDSS photo- 
metric survey catalogue is currently available for target se- 
lection. We use the published SDSS luminosity function to 
estimate the number of targets. In order to accurately es- 
timate the measurement error for the given sample, it is 
necessary to know the bias of the sample. We use the halo 
model to estimate the bias of the samples. 

It should be noted t hat in the methods employed here 
(|Seo fc EisensteinI |2003| . I2007D . the treatment of the non- 
linear effect and redshift distortion is very crude. Further- 
more, the assumption of scale inde pendent clustering bias 
may b e overly idealistic. Recently, ICrocce fc Scoccimarrol 
(' 20071 ) studied the BAO in both the power spectrum 
and correlation functio n by using the renormalized per- 
turbation theory(RPT, ICrocce fc Scoccimarrcl (|2006al lbl)). 
They showed that mode-coupling typically leads to percent- 
level shifts i n the acoustic peak of correlation function. 
ISmith et al.l (|2008l ) have also shown with both numeri- 
cal simulation and an alytical calculation that the peak is 
moved by these effects. lAngulo et al.l (|2008h proposed a new 
method to extract the unbiased estimation of sound hori- 
zon scale based on a large N-body simulation and semi- 
analytical galaxy formation model. ISeo et al~l (|2008l ) also 
presented high signal-to-noise ratio measurements of acous- 
tic scale from large volume simulations, and obtained a ro- 
bust measurements of acoustic scale with scatter close to 
that predicted by the methods utilized in this paper. There 
have also been works o n scale dependent halo and galaxy 
bias (jSmith et al.ll2007l ). In analysing the real data to ex- 
tract cosmological information, a more refined treatment is 
required to account for all these effects. 

We review the methods for error forecasting in § 2. 
§ 2.1 is devo ted to the Fisher matrix met hod, which was 
developed by ISeo fc EisensteinI (|2003l . l2007h for BAO mea- 
surements; § 2.2 is devoted to the MCMC method; in § 2.3 
we discuss how to estimate the bias of the sample for the 
planned surveys. After a brief description of the data that 
we assume LAMOST could collect, we present our forecast 
on the BAO measurement precision in § 3, and conclude 
in § 4. In the appendices, we describe the characteristics of 
the LAMOST telescope (Appendix A), and the details of our 
survey design (Appendix B), including the generic main sur- 
vey (Appendix Bl), the LRG survey (Appendix B2), and the 
QSO survey (Appendix B3). We estimate the density and 
bias of each sample, and also discuss the fiber time required 
for completing each survey. 

Throughout the pa per, we adopt the fl at ACDM, with 
the WMAP three year (ISpergel et al.ll2007^ best fltted pa- 
rameter values ( Q.bh^ = 0.0223, Q.mh^ = 0.128, Q.k = 
0, ft = 0.73, Wo = -1, wa = 0, = 0.958, As = 
2.3 X 10"'', T = 0.092 ) as our fiducial model. 



2 METHODS OF FORECAST 
2.1 Fisher Information Matrix 

Before recombination, the presence of large number of free 
electrons ensures that the photons and baryon plasma are 
tightly coupled. This provides a resilient force against any 
motion induced by gravity. Acoustic waves are generated 
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from the primordial fluctuations. When the Universe is suf- 
ficiently cooled for the photons and baryons to decouple, the 
oscillation ceases and the waves are imp rinted on the mat- 
ter and radiation density power spectra (IHu fc White 19961: 
Eisenstein et al.lll997l : lEisenstein fc Hu|[l998l : iMeiksin et. all 
19991 ). with the (comoving) characteristic scale (the BAO 
scale) determined by the sound horizon at the last scatter- 
ing surface: 



(1 + z)dt ■- 



Cs dz 
HIT)' 



(1) 



where Ca is the sound speed. For a given set of cosmologi- 
cal parameters, the absolute scale of sound horizon can be 
calculated. 

When making measurement of the matter power spec- 
trum in a galaxy survey, the distances along and perpen- 
dicular to the line of sight are measured from redshift and 
angular s eparation respectively. T he comoving distances are 
given by (ISeo fc Eisensteinll2003l ) 



{l + z)DAiz)Ae, 
cAz 



H{z) 



(2) 



The BAO scale provides a standard ruler to measure the 
angular diameter distance Da{z) and the Hubble expansion 
rate H{z). In a model with dark energy equation of state 
w(z), these are given by 

1/2 

n„(l- 



H{z) 



Ho 
and 

DAiz) 



1+ 
Jo 1 



(3) 



1 + 2 



dz 

mz)- 



(4) 



In this paper, we consider both a constant equation of state 
w for the dark energy, and a redshift-dependent one param- 
eterized in the form of 



io(z) = Wo + Wa[l — a.{z)] 



■ wo 



(5) 



l + z 

The statistical error in LSS measurements includes sam- 
pling variance due to the finite volume of the survey, as well 
as shot noise. In Fourier space, the statistical error is the 
summation of these two effects l|Feldman. Kaiser fcPeacockl 
ll994l : lTegmar3ll997l ): 



ap_ I 1 P{k) + 1/n ^ 

P ^ V Vsk^Ak P{k) ^ 



Vcaik)k^Ak' 



(6) 



where Vs is the survey volume and Vctt is the effective volume 
defined as 

2 



n(r-)P(fc,^) 



n{f)P{k,^l) + l 



dr. 



(7) 



Here jj, is the cosine of the angle between direction of k and 
the line of sight, and n{f) is the comoving number density 
of galaxies. For constant n. 



Vesik) 



nP{k) 



nP{k) + 1 



v.. 



(8) 



The Fisher information matrix for cosmologica l parame- 
ters derived from an LSS mea surement is given by l|Tegmarkl 
119971 : ISeo fc Eisensteinllioosl ) 



91nP(fe) ainP(fc),, dk 
Voff (fc) 



2(27r)3' 



(9) 



where pi denotes the cosmological parameters in the theory. 
In reality, the galaxy power spectrum is observed in red- 
shift space. The conversion from redshifts to physical scales 
depe nds on cosmological parameters (|Seo fc EisensteirJ 
[2003), 



-Pobs(^rcf_L 



Da{z 



,H(z) 



1+/3 



DAizyH,,t{z) 



xb^G{zfPn,oik)+Pa 



kl+kf. 



(10) 



where the subscript "ref" denotes quantities calculated in 
the reference cosmology, b is the bias factor for the galaxy 
sample, G{z) is the growth factor, and f3 is the redshift dis- 
tortion factor, k± and fcy are the k components perpendic- 
ular and parallel to the line of sight, respectively, and Pshot 
is the shot noise contribution. 

The measurement error of the power spectrum depends 
on the amplitude of the power spectrum, which is given by 
Pobs(fc) = b^Pun{k) in the k range relevant for BAO mea- 
surement. Thus, the uncertainty depends sensitively on the 
value of the bias b of the sample. We present our method of 
estimating the bias parameter in § 2.3. 

The measurements are also subject to some systematic 
errors, such as the effect of redshift-distortion due to peculiar 
velocity, non-linear evolution of the power spectrum, scale- 
dependent bias, and finite spectrograph resolution. 

Non-linear evolution of the density fluctuation en- 
hances small scale power, which corresponds to a smear- 
ing of the acoustic signature (Eisenstein, Scofc White 200^; 
[Seo fc Eisensteinll20(yi . This can be estimated by convolv- 
ing a Gaussian displacement field with the linear correlation 
function, 



Pnonl(fc, fl) = Plin(fc, m) SXp 



(11) 



where 



EoG, 



SoG(1 + /), 



with So = 12.4/i"^Mpc for as = 0.9, and / = dlnG/dlna 
is the derivative of the growth factor. The observed power 
spectrum is smeared below the resolution of the spectro- 
scopic measurement. We model such smearing as a Gaus- 
sian, 



P(fe, ^) = Pobs(fc, M) exp (-fey (Tr) 



(12) 



where ar — caz/H{z), and = {1 + z)ao, ctq = 1/P, with 
R the resolution of the spectroscopic measurement. We con- 
sider R — 1000 and R = 2000 respectively (for more details 
see Appendix A). The finite spectra resolution has the effect 
of smearing small scale powers. We set fcmax = 0.5/i/Mpc. 
The modes at very large scales contribute very little to the 
integral in equation (9), so we set kmin = 0. 

The parameters involved include cosmological parame- 
ters and survey parameters. To define a cosmological model, 
the following set of parameters are included: the Hubble 
constant h, baryon density Qb, matter density Qm, dark en- 
ergy density ilx, dark energy equation of state w, spectrum 
normalization In^s, spectral index Us, reionization optical 
depth T. In addition, for each redshift bin Zi of the survey, we 
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have parameters \nH{zi), In DA{zi), growth function G(zi), 
hnear redshift distortion \nl3{zi), and shot noise Pshot{zi). 

To obtain useful constraints on cosmological parame- 
ters, it is necessary to break the degeneracy by combining 
the BAO data with data obtained from some other cosmo- 
logical observations, e.g., CMB and/or Type la supernovae 
(SNIa). Consider first the combination of BAO and CMB, 
the total Fisher matrix is then given by 



(13) 



where F^^ {z„) is the Fisher matrix derived from the nth 
redshift bin of the large scale structure redshift survey. The 
CMB Fisher matrix is expressed as 



EE 



(C0v),^y- 



(14) 



where Cxi is the Z*** multipole for X — T (tem- 
perature correlation), E (E mode polarization corre- 
lation), B {B mode polarization correlation), and C 
(temperature-polariza ti on c ro ss-correlation), respe c tivel; 



temperature-polariza ti on c ro ss-correlation), respe c tively 
Seliak fc Zaldarriagal 1 19961: Kamionkowski et al.l 1 19961 : 



ISeliakHl997l :1 Eisenstein. Hu fc Tegmark||l999i r The elements 
of the covariance matrix Cov; between various power spectra 
are 

2 



(Cov)( JT 
(Cov);,_B£; 
(Gov);, SB 
(Cov)i, CC = 



[21 + l)/sky 
2 

(2/-H)/,ky 

2 

(21 + l)/,ky 

2 



{CEi+i^p'Bfy, 

{CBl+LOp'Bfy, 



{21 + l)/,ky 



{CIi + {Cti+uot^B^^) 



(Gov);, TC 
(Cov);,£;c 



(Cov);,TB 
2 



x{Cei+i^p^BI 
_ 2 

^ {21 + l)/,ky 



Ccu 



(15) 
(16) 
(17) 

(18) 
(19) 



{21 + l)/,ky 

2 



Cci{Cti+lot^B-^), (20) 
Cci{Cei+lop^B{^), (21) 



(2/ + l)/,ky 

(Cov),,TB = (CovObb = (CovOcs = 0, (22) 

where ujt,ojp are the inverse square of the detector noise 
level on a steradian patch for temperature and polarization, 
respectively. Bf = exp[— i(i + l)0bGain/8 In 2] is the beam 
function, ©beam is the full-width, half-maximum (FWHM) 
of the beam in radians. We shall assume that the CMB data 
would come from Planclfl, which is scheduled to be launched 
in 2009. By the time that the data of any of the surveys 
considered here is collected, Planck should have been in op- 
eration for at least three years. Therefore, we assume that 
CMB measurement errors correspond to Planck three-year 

observation (see Table [l]). 

In summary, similar to lSeo fc EisensteinI l|2003h , we con- 
sider the parameter set, 

{Pi,CMB} = {^mh? ,^bh^ ,^m,ns,\nAs,T,DA,CMB}, (23) 



Table 1. Assumed Planck Characteristic with three year integra- 
tion. 



Center Frequency (GHz) 


100 


143 


217 


Angular Resolution (arcmin) 


10.0 


7.1 


5.0 


At per pixel (pK) 


3.9 


3.5 


7.6 


Ap per pixel (jiK) 


6.3 


6.6 


15.4 


sky fraction /^ky 


0.65 


0.65 


0.65 



and the following additional parameters for LSS: 

{Pi,LSs} = {\nDA,iMHtMGiMPi,Pshot,i}- 

We then marginalize the nuisance parameters 
(In Gi, ln/3i, Ps/iof_i) by selecting the submatrix of F~.^ 
with appropriate column and rows. Finally, the Fisher 
matrix of dark energy parameters wo and Wa was ob- 
tained by converting from the distance parameter space 
{Pm} = {Dai, Hi} to the dark Energy parameter space 
{g,} = {wo, Wa, i^x}- 



pDB ^ dpm pdts dpn 

dq, ™" dqj 



(24) 



2.2 Monte Carlo Simulation 

Another method of making forecast is to use Monte Carlo 
to analyse a synthetic data set (e.g., power spectrum) 
with the expected error for the exp eriment in consideration 
(|Perotto et al.l l2006l : iLi et all 120081 ). This method is more 
time consuming than the Fisher matrix method, but it does 
not assume the likelihood to be nearly Gaussian. 

In order to deal with dynamical dark energy model, 
we modified the Boltzmann code camb togethe r with 
Monte Carlo program CosmoM(|f| (|Lewis fc Bridle! |2002| ) 
to search through the cosmological parameter space 
{Qth^ , i^ch^ , d, T, Wo, Wa,ns, As, {fv, rir) }, where 6 is the an- 
gle extended by the sound horizon at recombination and 
serves as an independent variable in lieu of the Hubble con- 
stant here. We generate the matter power spectrum and the 
CMB temperature and polarization angular power spectra. 
Again we assume Planck three-year ( Table. [1]) data that is 
generated by the program FutureCMBQ (|Perotto et al.ll2006l '). 
We do not include lensing effect here, since it has no signifi- 
cant impact on our results but requires more computer time. 

For the LSS part, we assume that the uncorrelated band 
power is a Gaussian realization with band error {ap/P)^ ~ 
l/(fcVcffAfc). The bin size of k is chosen to be similar to the 
SDSS-II LRG real data at low redshifts. The mock data in- 
clude larger range of k for higher redshifts to account for the 
fact that the linear approximation is valid at smaller scales 
then. The bias of the sample is also treated as a parameter 
in the MGMC. 

The overall constraint on dark energy could be im- 
proved significantly by adding data from complementary ob- 
servations. As we are interested in estimating what could be 
learned from "clean" signals such as BAO, and also what 



''' \http://www.rssd. esa.int/index.php?project=planck\ 



http:/ / cosmologist. info/ cosmomc/\ 

http://lappweb.m2p3.fr/r^perotto/FUTURCMB/home.htmi\ 
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could be achieved with aU types of experiments, we con- 
sider the constraints from two cases: (1) only the CMB and 
BAO data; (2) CMB, BAG and SNIa data. In the second 
case, we assume that a large numbe r of SNIa (1200 in total 
at < 2 < 1) collected by SNLS (|Astier et all |2006| ) are 
available. 



2.3 Estimation of the Bias Factor 

To make an estimate of the measurement with either the 
Fisher matrix or the MCMC method, one needs to know 
the bias of the sample, which determines the amplitude of 
the power spectrum. With the real data, t his can be mea- 
sured for the assumed cosmological model llCoil et al ]l2005l: 
Ross et al.ll2006l: iBlake. CoUister fc Lahavll2007l : IRoss et al.l 



Table 2. Summary of surveys. 



20071 : IPadmanabhan et al.ll2008l ). but in planning the survey 



one needs some way to estimate the value of the bias factor 
Here we use the halo mod el to make such an estimate. 
In the halo model Jjing et al. I 1 19981 : ISeliakI 



Peacockl I2OO0I: IScoccimarro et al.l I2001I : IShethI 
Coorav fc Shethll2002D r the number density of galaxies 
luminous than a threshold luminosity Lf 



dn 

dm — (N{m, > Lt)), 
am 



and the large scale galaxy bias factor is 



dn 

dm- — {N(m,> Lt))h(m) 
dm 



200c 



2001 



(25) 



(26) 



In the above expressions, dn/dm is the halo mass dis- 
tribution function, and b{m) is t he halo bias factor, fo r 
which we use the fitting formulae in ISheth fc Tormenl(|l999t ). 
The function {N{m, > L)) represents the halo occupation 
distribution function (HOD) for galaxies above luminos- 
ity L in halos of mass m, which can be parameterized as 
a sum of contributions from central and satellite galaxies 
JZehavi et al.ll2005bl : IZheng et al.]|2005l : iBerhnd fc Weinberg! 
I2OO2D . In our calculation of the galaxy number densities and 
large-scale bias factors, we adopt HOD parameters deter- 
mined from existing galaxy clustering data (see Appendix 
B). 

For t he large-scale bia s facto r of QSOs, we follow the 
model in IWvithe fc Loebl (|2003l ') (also see iMaruUi et aU 
I2OO6I ). In this model, the QSO activity is triggered by halo- 
halo major mergers. The black hole powering the QSO is 
assumed to be a fraction of the mass of the host halo. The 
QSO shines at the Eddington luminosity after the merger 
with a lifetime equal to the dynamical time of a galactic 
disk, with a universal light curve. The model gives a deter- 
ministic relation between the _B-band luminosity Lb and the 
halo mass m at a given redshift z and predicts the B-band 
luminosity function ^'(Ls,z). The bias factor for quasars 
more luminous than LB,min is given by 



b{z) 



b{LB{m,z))9{LB,z)dLB 



11 



^{LB,z)dLB 



(27) 



The model is not necessarily realistic, but for our purpose 
in this paper, we find that the uncertainty in the estimated 
bias from this simple model does not affect our results too 
much (see § 3.2). 



Sample 


Magnitude 
Limit 


Surface Density 
(deg-2) 


Targets 
(106) 


Fiber Hours 
(106) 


MAINl 


r < 18.8 


330 


2.6 


0.88 


MAIN2 


r < 19.8 


1050 


8.4 


14 


LRG 


idev < 20 


205 


1.5 


1.9 


QSOl 


g < 20.5 


30 


0.24 


0.06 


QS02 


5 < 21 


45 


0.36 


0.18 


QSOS 


g < 21.65 


72 


0.57 


0.855 



3 RESULTS 

We shall consider three surveys for the LAMOST telescope: 
(1) a magnitude-limited general survey of galaxies of all 
types which we shall call it the main survey; (2) an LRG 
survey; (3) a magnitude-limited low redshift quasar survey 
{z < 2.1). The characteristics of the LAMOST telescope and 
the details of the sample construction and survey design are 
discussed in the Appendices. 

Table [2] gives a summary of these surveys, including 
the magnitude limit, the surface density, the total number 
of targets, and the required fiber hours, i.e., the number of 
fibers allocated for the survey multiplied by the observing 
hours. For the main and QSO surveys, several different mag- 
nitude limits are considered. A general galaxy survey with 
r < 18.8 (denoted as MAINl) may be completed within one 
year, while a deeper one with r < 19.8 (denoted as MAIN2) 
may serve as the ultimate goal of the LAMOST telescope 
operation. For the QSOs, three magnitude limits are con- 
sidered, g < 20.5 (QSOl), g < 21.0 (QSG2), and g < 21.65 
(QSG3), aU between z = 0.4 and z = 2.1. For the LRG sur- 
vey, we consider a sample tha t is similar to the SDSS MegaZ 
sample (|Gollister et. al.ll2006l ). Based on the observationally 
inferred luminosity function and HOD parameters, we esti- 
mate the comoving density n[z) and the clustering bias b 
to determine the effective volume defined in Eq. ((T)). Since 
the angular density of LAMOST fibers is 200 deg"^, these 
targets may not be observed in a single run. We discuss in 
more detail the fiber allocation strategy in the Appendices. 

The redshift distribution of galaxies is shown in Fig. [1] 
(the comoving densities are shown in Fig. IB3|) . We see that 
the redshift distribution of the MAINl and MAIN2 samples 
are broader than the SDSS main sample. The minimum red- 
shift of the LAMOST LRGs considered here is 0.4. It may 
be possible to have some more low redshift LRGs, but the 
present one is based on a sample which has already been 
selected (MegaZ). 

The redshift distribution of quasars is shown in Fig. [2] 
(comoving density is shown in Fig. IB4|I . Compared with the 
galaxies, the number of quasars is much smaller. The drop 
at low redshift is a combined effect of the existence of lower 
absolute magnitude limit (below it an object would be classi- 
fied as an AGN rather than a quasar) as well as the decrease 
of quasar activity. 



3.1 Effective Volume and Distance Scale 

In Fig. Owe compare the effective volume of different LAM- 
OST survey samples considered in this paper (LAMOST 
MAIN, LRG and QSO) and the SDSS samples (SDSS main 
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Figure 1. The Redshift distribution of diflferent galaxy samples, 
assuming a sky coverage of 8000 deg^ 
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Figure 2. The redshift distribution of the QSO samples. 

and LRg13). While at small k (large scales) the effective 
volumes of the galaxy samples are almost constant, at small 
scales their effective volumes decline. The decline occurs at 
a scale where shot noise starts to dominate over the power 
spectrum (Eq. |8]). So the effective volume of a sparsely dis- 
tributed sample drops at a scale larger than that of a high 
density sample. The finite spectrograph resolution smears 
the power spectrum at small scales and further reduces the 
signal. We plot the effective volume for both R = 1000 (solid 

For the SDSS LRG, we assume th at the survey volume i s 
0.8(Gpc//i)3 for 8,000 deg^, whereas [Seo fc EisensteinI ||2003| ') 
assumed a volume of l(Gpc/h)'' for 10000 deg^. The density 
is assumed to be 10~*(/i/Mpc)^, an d for the bias we ta ke 
the measured va l ue (1 .9, as given in iTegmark et al.l ll2006f) ). 
iBeo fc EisensteinI ||2003|') originally assumed a bias of 2.13. The 
actual s urvey data may ext e nd to h igher redshift, see lEisensteinl 
ll2005d) : [CabrefcGaztanaeal ll2008al lbh. 
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Figure 3. Effective volume of different surveys. Solid curves for 
resolution R = 1000, dotted curves for R = 2000. 

curve) and R — 2000 (dotted curve) in the figure, and we see 
that the decline in effective volume starts at smaller scale 
for the higher resolution setting. 

As shown in the figure, the MAINl survey has an ef- 
fective volume of 1 (Gpc//i)^, about 3.3 times that of the 
SDSS main survey on large scales - it is even larger than 
the SDSS LRG survey as assumed in this paper (though the 
actual SDSS LRG survey reaches a greater effective volume, 
see footnote 10). If the survey of this sample is completed 
within one year as planned, a second round survey could 
go one magnitude deeper and obtain the MAIN2 sample, 
which has the large effective volume of 3 (Gpc/h)^. However, 
as discussed in Appendix B, collecting this sample would re- 
quire much more observing fiber time, since both the number 
of targets and the integration time required for each target 
increase. 

If an LRG sample similar to that of the MegaZ catalog 
(but extended to 8, 000 deg^ of the sky) is surveyed, the ef- 
fective volume at large scales is 3.1 (Gpc//i)^, 4 times larger 
than the SDSS LRG sample. While the LAMOST LRG sam- 
ple has an effective volume comparable to that of LAMOST 
MAIN2 sample, the LRG survey would take much less time. 

The shape of the QSO effective volume is different from 
those of galaxies. At large scales, while the effective volume 
of a galaxy sample becomes flat, that of the QSO sample 
increases with k. This is because the density of QSO is very 
low and the shot noise always dominates. As a result, the 
effective volume traces the shape of the power spectrum (c.f. 
Eq. m with nP(fc) < 1). The effective volume of the QSO 
sample can reach a peak which is nearly 3 times larger than 
that of the LRG sample at 0.02/i/Mpc, thanks to the large 
survey volume. However, it decreases quickly towards small 
scales, and becomes comparable to the SDSS main sample at 
k ~ O.lft/Mpc. Since the decline happens at scales greater 
than the BAO scale, the constraining power of the QSO 
sample on cosmology is relatively weak, although we still 
expect to detect the BAO signature with the QSO sample. 
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Figure 4. The expected fractional error of tlie galaxy and quasar 
power spectra. 
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Figure 5. The distance scale Dji{z) and H{z) for the fiducial 
model and the expected measurement errors. The first LRG point 
is omitted due to coincidence with other data points. 



In the above we have discussed the total effective vol- 
ume of the sample. In practice, the survey volume will be 
divided into several redshift bins. The sample density is ap- 
proximately a constant within each bin, and the cosmolog- 
ical distance scale is measured at the center of each bin. 
The measurement errors on the power spectrum and dis- 
tance scale depend on the size of the redshift bins. Larger 
bins lead to smaller sample variance but fewer data points 
on distance measurements, and vice versa. Our tests show 
that the eventual errors on dark energy equation of state or 
other cosmological parameters are not sensitive to how we 
make the division of the bins, as long as the evolution in 
density and bias within each bin is small. 

In Table |3] we show our division of redshift bins for 
different samples. For the main samples, our bin size is 
Az = 0.2, and the MAINl and MAIN2 samples have two and 
three bins, respectively. For the LRGs, the first bin contains 
the SDSS LRG sample, and the other two bins are similar 
to those in the MegaZ catalogue. The three QSO samples 
have different turnover redshifts (determined by the appar- 
ent magnitude limit as well as the absolute magnitude limit) 
and they are divided into slightly different redshift bins. The 
bias of each sample is given in Table [B2l 

We plot the forecasted matter power spectrum measure- 
ment errors derived from these samples at different redshifts 
in Fig. m For comparison, we also show the expected errors 
for the SDSS LRG sample. These errors are calculated as- 
suming that the whole survey volume is included. For the 
LAMOST galaxy samples, the fractional error is at a level of 
1% around k = 0.1/i/Mpc. The errors for the QSO samples 
are large on small scales due to their low space densities, but 
on large scales {k ~ 0.01/i/Mpc) QSO samples yield smaller 
errors than galaxy sample because of their large survey vol- 
umes. 

In Fig. Owe plot the distance scales {Da, H) and their 
measurement errors {odai ^h) as a function of redshift for 
the fiducial model. The top and bottom panels show the 
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Figure 6. Constraints (Icr) derived from combing different data 
sets for a dark energy model with constant equation of state ui, 
all assuming the Planck 3-year prior. 



errors for H{z) and Da{z), respectively. The measurement 
errors for the galaxy samples are at the level of a few percent. 
The QSO samples have larger measurement errors, but they 
can extend to higher redshifts. 



3.2 Constraints on Cosmology — Results from the 
Fisher Matrix Method 

Using the projection method described in §2.1, we derive the 
error ellipse on the dark energy parameters after marginal- 
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Table 3. Redshift bins of different surveys. 
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Figure 7. Constraints on dark energy parameters wq and Wa de- 
rived from combining different data sets, all assuming the Planck 
3-year prior. 



izing over other parameters. In Fig. |6l the constraints on 
a dark energy model of constant equation of state w are 
plotted for the MAINl sample and the combination of the 
MAINl, LRG and QS03 sample. For comparison, we also 
plot the expected error ellipse from the SDSS LRG sample, 
which has a good constraint on the dark energy density pa- 
rameter Q.X and a poor constraint on w. With the LAMOST 
samples, even for the MAINl only sample, there is a signifi- 
cant improvement in constraining Q.x and w. By combining 
the MAINl sample, the LRG sample, and the QSO sam- 
ple, the allowed range in w is halved with respect to that 
from the SDSS LRG, reaching a fractional uncertainty of 
about 6.4%. We did not plot the constraint produced by the 
MAIN2 sample in the above figure, but it is almost the same 
as the combination of MAINl, LRG and QS03 sample. 

Next we consider the case of evolving equation of state 
w, which is common for dynamical dark energy models. In 
Fig. [T] we show the constraints on dark energy equation of 
state parameters from different combinations of the MAIN, 
the LRG, and the QSO samples. The solid contour shows the 
dark energy constraints that the LAMOST can achieve with 
all the samples. There are significant improvements over pre- 
vious surveys such as the SDSS. However, the constraints 
on evolving w are still weak, which is in agreement with 
the assessement of t he Dark Energy Task Force (DETF) 
lAlbrecht et al.ir2006l ) that near term projects could not yet 



Table 4. uiq — Wa FoM for different samples with R = 1000 and 
R = 2000. The fourth column shows the fractional increasing of 
FoM by raising the resolution. 



Sample/FoM 


R=1000 


R=2000 


+% 


MAINl 


5.72 


5.79 


1.3% 


MAIN2 


10.83 


11.04 


1.9 % 


LRG 
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8.17 


2.9% 


QSOl 


0.554 


0.638 


15.2% 


QS02 


0.903 


1.049 


16.2 % 
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1.425 
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17.4 % 





Figure 8. The Icr error ellipse for too aid Wa for different sam- 
ples: the main samples (top left), LRGs (top right), QSOs (bot- 
tom left), and comparison of main and LRGs (bottom right), all 
assuming the Planck 3-year prior. 



provide high precision measurement on uiq and Wa- By com- 
bining the BAO measurements with future Type la super- 
novae results, the constraints could be greatly tightened, as 
shown by the dotted contour in the plot. 

In the above we have considered the constraints de- 
rived by combining samples of different LAMOST surveys. 
How much information does each of the samples provide? 
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To address this problem we plot the constraints obtained 
with each of the samples in Fig. [8] We can see that the 
galaxy samples provide stronger constraints than the quasar 
samples, the constraining power of the LRG is between the 
MAINl and MAIN2 samples. 

We can also define a Figure of Merit (FoM) 



FoM = 



1 



(28) 



•\/det(Cov(wo,ma)) 

which is proportional to the inverse of the contour aref'l. 
It is 2.01 for the SDSS LRG sample, 5.72 for the MAINl 
sample, 7.94 for the LAMOST LRG sample, 10.83 for the 
MAIN2 sample. The QSO samples all have very large survey 
volumes, but due to their low densities, their FoMs are small: 
FoMqsoi = 0.55, FoMqso2 = 0.90, and FoMqsos = 1-43. 
The error ellipses in Fig. [S] reflect this trend: the MAIN2 
sample provides the most stringent constraint and the LRG 
observation could also achieve a similar FoM with much less 
fiber time. When combining these samples, we reach a FoM 
of 20.88. Note that since the MAIN2 sample overlaps the 
LRG sample, we only include the first two bins of MAIN2 
in the combination. 

The estimate on the FoMs for galaxy samples is not 
very sensitive to the adopted value of the bias factors, as 
long as the constraining power comes mainly from the range 
of the power spectrum that is not dominated by shot noise 
(Eq. |8}. We test this by adopting observationally inferred 
bias factor, instead of that from the HOD method. If we 
ignore the difi'ere nce in the sample definition and use the bias 
facto rs b = 2.35 ||Ross et al.ll2007l ') and b = 1.66 (|Ross et all 
bOOel ) for the LRG sample, we find that the FoM of the 
LRG sample changes by 1.8% and -6.8%, respectively. For 
the QS03 sample, we find that th e change in th e FoM is 
6.9% if the bias factor measured by iPorciani et al.l (|2004l ) IS 
used in our analysis. 

An important question in the design of the surveys is 
whether to choose R = 1000 or R = 2000. With the higher 
resolution setting, the redshift could in principle be mea- 
sured to higher precision. However, in practice, for a fixed 
observation time, a higher signal to noise ratio could be 
achieved with the lower resolution setting. Conversely, if 
one requires a fixed signal to noise ratio, more time would 
be required to complete the survey for the higher resolu- 
tion setting. For some applications, particularly those on 
the physics of the galaxies, useful information could be ob- 
tained only with the higher resolution setting. For dark en- 
ergy constraint, the difference between the two choice is 
largely quantitative. From Fig. [3] we see that there is al- 
most no difference on large scales between the two resolu- 
tion settings. On small scales, however, the effective volume 
is larger for the higher resolution setting, which shifts the 
smearing of small scale power to smaller scales. From the 
figure, it seems that for galaxy samples the improvement is 
larger than the QSO samples. However, for the QSO the im- 
provement from the higher resolution occurs at the critical 
scale of about 0.05 ~ 0.1/i/Mpc, which is exactly where the 



Our definition follows that o f lWand 1 I2OO8I) . In tlie Fisher ma- 
trix formalism, it is related to the DETF FoM, which is the in- 
verse of the 2a contour area llAlbrecht et ahlboOSh as FoM^ang = 
6.177r FoMdetf- 
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Figure 9. The probability distribution function for some cosmo- 
logical parameters constrained by combining the CMB and BAO 
data. The solid and dashed curves correspond LAMOST LRG 
(MegaZ) sample and the SDSS LRG sample. The 68.3% distribu- 
tion is marked with two verical lines. 



baryon oscillation features occur. The figure of merit of dark 
energy measurement is given in Table |4] for the two settings. 
For the QSOs, the FoM with R = 2000 is about 15% bigger 
than R = 1000. For the galaxies including both the MAIN 
and LRG samples, the difference is only a few percent, since 
the improvement is made largely in the non-linear regime. 
Therefore, from the perspective of dark energy constraints, 
R — 1000 is sufficient and preferred. 



3.3 Constraints on Cosmology 
MCMC Method 



Results from the 



We have also performed Monte Carlo simulations for a few 
cases. In particular we consider the SDSS LRG, which repre- 
sents our currently available data, and the LAMOST LRG, 
which represents the precision achievable after two or three 
years of LAMOST survey. We derive constraints for two 
cases: (1) only the CMB and BAO data are used; (2) CMB, 
BAO and SNIa data are used. In the first case, only the 
"clean" methods (i.e., the CMB and BAO) that are based 
on well understood physics are used. In the second case, we 
include a sample of 1200 supernovae expected from SNLS. 

In Fig. |9]we show the probability distribution function 
(PDF) for each of the interested cosmological parameters ob- 
tained from the first case after marginalizing over the rest 
of the parameters. The solid (dashed) curves are from CMB 
combined with the LAMOST LRG (SDSS-II LRG) sample. 
Except for those of wq and Wa, the PDF is well approximated 
by a Gaussian distribution. There is no significant improv- 
ment in the constraints on cosmological parameters other 
than Wo and Wa by the addition of the LAMOST data, as the 
constraints from the CMB data are already very stringent. 
Compared with the SDSS-II LRG sample, the LAMOST 
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Figure 10. The expected constraints on dark energy equation 
of state parameters by combining the CMB data and the BAO 
measurements with the LRG sample. The BAO measurements in 
the top and bottom panels are assumed to come from the SDSS-II 
LRG and the LAMOST LRG samples, respectively. Contours in 
each panel correspond to 68.3% and 95% confidence levels, and 
the shading shows the mean likelihoods. 



LRG sample helps to significantly tighten the constraints 
on the dark energy parameters. 

In Fig. [TO] we show the constraints in the Wo-Wa plane. 
The contours correspond to 68.3% and 95% confidence lev- 
els, and the shading is the mean likelihoods of samples 
l|Lewis fc Bridlell2002( ). With the LAMOST LRGs, the con- 
tours are much smaller than those with the SDSS-II LRGs. 
For our choice of non-linear cut-off scale, the size of the con- 
tour generally agrees with the Fisher matrix estimation, but 
not exactly. In particular, for SDSS-II LRGs, the long tail 
towards the lower right corner reflects the non-Gaussinity of 



Figure 11. Top: Constraints on uiq and Wa by combining the 
CMB (Planck) and the SNIa (SNLS) data. Bottom: Constraints 
on Wo and Wa by combining the CMB (Planck), SN la (SNLS) 
and also the BAO (LAMOST LRG) data. Contours correspond 
to 68.3% and 95% confidence levels, and the shading shows the 
mean likelihoods. 



the likelihood distribution. With only the CMB and BAO 
data measured by LAMOST, the constraint on dark energy 
is not very stringent, and many different dark energy model 
would have easily passed the test. 

To improve the constraints on dark energy parameters, 
we can bring in the SNIa data. To separate the contributions 
from the SNIa and BAO data, we first show the constraints 
with CMB and SNIa data in the top panel of Fig. [11] Com- 
paring with the top panel in Fig. [10] (CMB-f BAO), we see 
that the SNIa data provide significant complementary infor- 
mation. The bottom panel of Fig. [TT] shows the constraints 
from all the three data sets (CMB-|-SNIa+BAO). Compared 
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with CMB+SNIa, adding BAO data from LAMOST LRG 
further tightens the constraints, increasing the FoM from 
33.07 to 42.83. Here, the FoM is defined as 6.17tt over the 
contour area of 95% confidence level. Such constraints may 
be able to distinguish some dynamical dark energy models. 



4 CONCLUSION 

In this paper we forecast the measurement precision of dark 
energy equation of state from the BAO measurement ex- 
pected from LAMOST surveys. We investigate three types 
of sample targets: magnitude limited surveys of all types of 
galaxies (main survey), a survey of color-selected luminous 
red galaxies (LRG), and surveys of quasars (QSO). For the 
main survey and the QSO survey, we consider a few different 
magnitude limits, all deeper than the SDSS spectroscopic 
survey. For the LRGs, we have used the MegaZ sample as 
our reference sample. The required fiber hours to complete 
each survey are also estimated (Tabl^. 

For each sample, we predict the redshift distribution, 
and use the halo model to estimate their clustering bias. We 
then calculate the effective volume of each survey as a func- 
tion of wavenumber k, the statistical errors in the measured 
matter power spectrum, and the resulting errors in distance 
scale measurements with the BAO used as standard rulers. 
The constraints on the dark energy equation of state pa- 
rameters are derived using two methods, the Fisher matrix 
method and the MCMC method. 

Our results yield useful insights on how to design the 
surveys. We find that the MAINl sample, which is one mag- 
nitude deeper than the SDSS main sample has an effective 
volume about three and half times larger than the SDSS 
main sample. Similarly, with the LAMOST LRGs, an ef- 
fective volume of about four times of the SDSS LRG can 
be obtained. This effective volume is almost as large as 
the much more time-consuming MAIN2 sample (two magni- 
tudes deeper than the SDSS main sample). The QSO sam- 
ples can also yield high effective volume at large scales, but 
it declines rapidly at small scales. The improvement on the 
measurement of dark energy equation of state depends not 
only on the large scale effective volume, but also on the ef- 
fective volume at relatively small scales. Our investigation 
shows that although the effective volume of the LAMOST 
LRG sample is comparable to the MAIN2 at large scales, 
the latter's figure of merit for dark energy measurement is 
greater. The QSO sample, on the other hand, generally has 
lower figure of merit. However, the number of targets and 
required observing fiber hours of the QSO samples are also 
much smaller than the galaxy samples, so it should not be 
regarded as a competing project of the galaxy surveys. The 
efficient design for using LAMOST survey to constrain dark 
energy parameters is to have a MAINl survey, an LRG sur- 
vey supplemented by a QSO survey. 

We also studied the impact of spectral resolution. The 
figure of merit improves by only about l%-3% for the galaxy 
samples when one switches the spectral resolution from 1000 
to 2000. For quasars the figure of merit improvement is 
about 15%. From the point of view of dark energy con- 
straints, it is unnecessary to achieve a high spectral reso- 
lution, although high resolution spectra are useful for other 
scientific objectives. 



Finally, we also used Monte Carlo simulations to make 
the same forecast to test the accracy of the Fisher matrix 
result. The figure of merits of the two estimate roughly 
agree, although the contours obtained from the Monte Carlo 
method are not exactly ellipses, but rather asymmetric. 

Overall, the LAMOST surveys could achieve an figure of 
merit (as defined bv lWand (120081 ) ) of 5-10 from each individ- 
ual galaxy samples. For a constant w dark energy model, the 
expected precision on w is about 6.9%. For time-dependent 
w models, the LAMOST survey would also help to signifi- 
cantly improve the constraints on the equation-of-state pa- 
rameters, although the expected errors on wo and Wa are 
not tight enough to pin down the nature of dark energy. 
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APPENDIX A: THE LAMOST TELESCOPE 

The LAMOST telescope (see e.g.. IChu|[l99i : IStondl200"i ) is 
a 4 meter telescopj^. its design features an fixed horizon- 
tal meridian refiecting Schmidt configuration, with an active 
Schimidt correcting plate to achieve both wide field of view 
(5° diameter, or 20deg'^) and large aperture. The focal plane 
accommodates 4000 optical fibers for spectrographic obser- 
vations. The optical axis is fixed in the meridian plane, tak- 
ing measurements of targets as they pass through the merid- 
ian (typically a target can be tracked for about 2 hours.) 
The telescope is located at the Xinglong station (117°34'30" 
East, 40° 23' 36" North) of the National Astronomical Obser- 
vatory of China, and can observe the part of sky with decli- 
nation -10° < 6 < 90°. The construction of the LAMOST 
telescope was completed by August 2008 and it is undergoing 



The main reflector of the LAMOST has a diameter of 6 meters. 
The mirrors of LAMOST are hexagonal, the efli'ective aperture for 
the whole telescope varies between 3.8 m to 4.5 m depending on 
the position of the pixel on the focal plane and the altitude of the 
target. Calibration for fibers on different part of the focal plane 
Is provided by observing standard stars at real time. 
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Table Al. Main Characteristics of LAMOST 

Average aperture 4 m 

o 

Field of view 5 

Focal plane (^1.75m 

Focal length 20 m 

Number of fibers 4000 

Spectral ranges 370-900nm 

Spectral resolution 1000 or 2000 

Minimal distance of adjacent targets: 55' 

Observable sky -10 ° < 5 < +90 ° 



engineering tests and calibrations. The surveys are expected 
to begin in 2010/2011 and last for at least 5 years. 

The LAMOST is to be equipped with both low reso- 
lution spectrographs and medium resolution spectrographs 
l|Zhu fc XulbOOol : IZhu et al.ll2004) , the latter is used mainly 
for bright targets such as stars. For extragalactic observa- 
tions using the low resolution spectrographs, spectra can be 
taken with two resolution settings: (1) using the full sized 
fiber, the resolution is 1000 (2) using half sized fiber, the 
resolution is 2000. The main features of the telescope is sum- 
marized in Table IXT] 

The LAMOST has a parallely c ontrolled fiber posi- 
tionin g system with 4000 optical fibers l|Xull2003l : IZhu et al] 
l2004h . each individually driven by two micro-stepping mo- 
tors and can move within a circle of 30 mm diameter (cor- 
responding to 340' on the sky) to aim at a target within the 
circle. The distance between two adjacent circle axes are 26 
mm, slightly smaller than the diameter of the circle, so that 
adjacent circles overlap. This allows full coverage of the fo- 
cal plane. Due to fiber collision, targets within 55 arcsec can 
not be observed at the same time. The average angular den- 
sity of the fibers is about 200 per square degree. The whole 
assembly is controlled by computer and can be adjusted in 
real time. Acquiring new target and adjust the fibers takes 
a few minutes typically. 

As the LAMOST telescope is dedicated to spectroscopy 
and does not have imaging capability, photometric catalogue 
compiled from other surveys would be used to select tar- 
gets. The LAMOST will carry out several different survey 
projects simultaneously. The input catalogue of these differ- 
ent surveys will be combined (with different assigned prior- 
ity) to produce t he general input catalogue of LAMOST. A 
software system l|Peng et al.ll2004r ) would then select targets 
of observation from this catalogue, taken into account the 
priority of the target, the observing time and target area, 
the lunar phase and weather condition etc, this allows more 
efficient use of the fibers and observing time. During each 
run, the selected targets are a mixture of different survey 
projects, so the appropriate measure of observing resource 
is the fiber time, i.e. the number of fibers allocated multi- 
plied by the observing time. 

The collecting area of LAMOST is about 2.2 times that 
of the SDSS telescope, which has an aperture of 2.7m. We 
would expect that for comparable exposure time, target of 
one magnitude deeper could be observed if the system effi- 
ciency and sky background are the same. The system effe- 
ciency, in turn, also depends on the positioning error of the 
fiber. The signal to noise ratio of an observation with differ- 
ent conditions and parameters has been estimated using ray 
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Figure Al. Expected exposure time of different bands and model 
profiles, assuming spectral resolution 1000, S/N=4 per spectral 
pixel, a sky background of 20.5 in the respective band (r or i), 
fiber positioning error 0.5", re = 3" for exponential profile and 
Te = 4" for de Vaucouleurs profile. 

tracing m ethod, and a calculator has been developed for es- 
timation (jXue 11200^ 1. For the present paper, we assume an 
average fiber positioning error of 0.5", and a sky background 
of 20.5 in r band We also assume that an average signal 
to noise ratio of 4 per spectral resolution pixel is achieved 
for both galaxies and point sources. This is about the min- 
imum for redshift measurement. Surveys with low signal to 
noise ratio allows fast survey of a large number of targets, 
however, in the real survey, a higher signal to noise ratio 
may be required for other scientific purposes. 

The required exposure time for achieving the desirable 
signal to noise ratio depends also on the size of the object 
(characterized by the half-light radii of the image 7-50 for 
extended source), the seeing of the sky, and the point spread 
function of the optical system. We characterize the angular 
size of the target by a parameter ■ For the de Vaucouleurs 
profile which represents elliptical galaxies, 

/dcv(r-) = /oexp[-7.67(r/re)'/*], 

and for the exponential profile which represents spiral galax- 
ies, 

/oxp(r-) = /oexp [-1.68r/re]. 

Analysis o f the half-light rad ii ij-^o) of SDSS galaxies 
show that (|Blanton et al.l 1200 ll ) about 80% galaxies have 
rso.petro ^ 3". For the exponential profile r5o,petro = 
I 50, true — Te, and for the de Vaucouleurs profile, rso, petro ~ 
0.77-50, true = 0.7re. Thus, for a simple estimate of the size of 
the targets, we assume that the average target has re = 3" 



^ Recent observations at Xinglong station show a sky back- 
ground between the 20.5 and 21.0 during moonless nights (H. 
Wu, pravite communication). 
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and 4" for exponential profile and de Vaucouleurs respec- 
tively. 

In Fig. lAll we plot the exposure time required for tar- 
gets of different sizes, bands and model profiles. As can be 
seen from the figure, for a given S/N, the integration time 
is proportional to for bright sources, where / is the 
flux, and to for the faint sources where the sky back- 
ground noise dominates. Compared with elliptical galaxies, 
the spiral typically requires more time for observation to 
reach the same magnitude, this comprises about 25% - 50% 
of our sample. Observing large samples of targets fainter 
than r — 21 would be impractical. 

The spectral resolution of the low-resolution grating is 
1000 when the full sized fiber is used, 2000 when half sized 
fiber is used. When observing a galaxy or quasar with mul- 
tiple lines, the redshift could be measured with accuracy 
better than the grating resolution if multiple spectral lines 
are present and the signal to noise ratio is high. However, 
the number of usable lines may be limited for many galaxies, 
and to save observing time the signal to noise ratio is per- 
haps not very high in large scale observations, so we assume 
the redshift measurement precision to be 1/1000 or 1/2000. 



APPENDIX B: DESIGN OF THE 
EXTRAGALACTIC SURVEYS 

Since the LAMOST does not have imaging capability by 
itself, the targets of the spectroscopic survey must be se- 
lected from photometric catalogues compiled from observa- 
tional data obtained by other instruments. At present, the 
SDSS photometric catalogue is the largest of such catalogue, 
which covers about 8000 deg'^ of northern sky up to r = 23, 
and only a small fraction of the objects on this catalogue 
has been targeted by the SDSS spectroscopic surveys. One 
may also consider to conduct a photometric survey for a 
part of the south galactic cap to supplement this. Regions 
of low galactic latitude are of less interest for extragalac- 
tic surveys due to the higher extinction. In the future more 
targets could be obtained from photometric surveys which 
are being planned, e.g., the Pan-Starrs surveys. However, as 
discussed in Appendix A, it is unlikely for LAMOST to ob- 
serve a large number of very faint objects within the limited 
observation time available. 

In the present work, we shall consider three surveys: (1) 
a magnitude limited general survey of galaxies of all types 
which we shall call the main survey; (2) an LRG survey; (3) 
a magnitude limited low redshift {z < 2.1) quasar survey. In 
Table [5] we list a summary of these surveys, including the 
surface density, total number of targets, and required fiber 
hours. 



Bl The main survey 

In the original feasibility study of the LAMOST, a typical 
number quoted is B=20.5 for an exposure time of 1 hour, 
this is comparable to the present estimate. The SDSS team 
has chosen their main sample using an r-band limit, as the r- 
band luminosity is a better tracer for mass. Here we assume 
the survey limit is given in r-band. The magnitude limit we 
quote here is just a strawman value and can be changed 
for the real survey. It is reasonable to consider a sample 



which is one or two magnitudes deeper than the SDSS main 
sample. The SDSS main sample has a magnitude limit of 
r=17.77, the exposu re time is typically about 45 minutes 
jStrauss et al.ll2003l ). We consider two samples: a magni- 
tude limited sample which is one magnitude deeper than 
the SDSS main sample, i.e. r < 18.8 (hereafter MAINl), 
and also a sample which is 2 magnitudes deeper than SDSS 
main, i.e. r < 19.8 (hereafter MAIN2). 

Given the above magnitude limit, we may estimate 
the angular density of targets using galaxies number count 
Yasuda et al.l 2001 ) or the observed lum inosity function 



Blanton et al.l 



20031 ). I Yasuda et all l|200ll ') have shown that 



the observed number of galaxies N(m) can be well fitted 
by assuming the Euclidean geometry for three-dimensional 
space, which gives the number of galaxies brighter than m 



N(m) oc 10° 



(Bl) 



Alterna tively, the luminosity function of SDSS galaxies was 
given in lBlanton et al.1 (120031 ) 

$(M,z) = 0.41nl0x$n0°-*'"+i)(*^*-*^' exp(-10°-^(^^*-"')(B2) 

and the redshift evolution of the luminosity is modelled by 

$*(z) = $*(0o)lO°-^^("-"«> (B3) 
M*{z) = M" {zo) - Q{z - zo) (B4) 

The absolute magnitude M is related to the apparent mag- 
nitude m and redshift z as 



M - 5 logio h = m~ DM{z) - K{z) 



(B5) 



where DM{z) is the distance modulus and K{z) the 
K-correction. The K-correction depends on the SED of 
the galaxy, we have adopted a simple approximation of 
K{z) ~ —0.21 -|- l.lz, which is approximately the me- 
di an K-correction for the 12 types of galaxy SEDs given 
m iBlanton et ai] l|2003h . This may not be very accurate for 
individual galaxies, but should suffice for the purpose of es- 
timating the angular density. 

We show in Fig. IBll the integrated angular density as a 
function of the magnitude limit in r band. The solid curve 
is obtained by usin g the galaxies number count given in 
lYasuda et al.l (|200lh . and the dashed curve i s obta ined by 
using the luminosity function of lBlanton et al.l (|2003l ). These 
data have to be extrapolated to faint end or high redshift 
to obtain the distribution function. We use the best-fitted 
r band evolut i on pa rameter Q = 1.62, P = 0.18 given in 
iBlanton et all (|2003l '). In this paper, we have used the num- 
ber count directly for the fibre time estimation, but for cos- 
mological constraints we used the luminosity function. There 
is only slight difference between these two. 

For the MAINl (r = 18.8) sample, the angular density 
is 330 deg"'^, about 1.65 times the average angular density 
of the fibers. Most targets can be observed in one to two 
observing runs (due to fiber collision, there could be a small 
fraction which is left unobserved). For MAIN2, the angular 
density is about 1050 deg"'^. These can be observed in 5-6 
observing runs. 

The total number of target is about 2.6 million for the 
MAINl survey, and 8.4 milhon for the MAIN2 survey. The 
total fiber time required for completing the survey is about 
0.88 million fiber hours for MAINl, and 14 million fiber 
hours for MAIN2. If we assume that all fibers are allocated 
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Figure Bl. The angular density of objects as a function of lim- 
iting magnitude in the general galaxy survey. The solid line cor- 
responds to the (apparent) magnitu de integration of the galaxies 
number count ijYasuda et al.l l200lh at SDSS commission stage. 
The dashed line is the result of (absolute) magnitude-spatial inte- 
gration of the luminosity function (Blanton et al. 2003), by using 
the best-fitted r band evolution parameter Q = 1.62, P = 0.18 
and a simple K correction model. 



Table Bl. HOD parameters for galaxies (reproduced from 
IZehavi et all ll2005d) ). 
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-21.0 
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14.09 


1.39 


-20.5 
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-20.0 


12.01 
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1.16 


-19.5 


11.76 


13.15 


1.13 


-19.0 


11.59 


12.94 


1.08 


-18.5 


11.44 


12.77 


1.01 


-18.0 


11.27 


12.57 


0.92 



for this survey, and for each night the observing time is 
6 hours, then MAINl observation can be completed in 37 
clear dark (moonless) nights, and MAIN2 in 584 clear dark 
nights. In reality, the automatic survey strategy system will 
probably assign targets of different surveys in the same sky 
area at the same time, and only 1/3-1/2 of the nights have 
good weather. Thus, we expect that the MAINl sample can 
be surveyed in the first year, while the MAIN2 sample in the 
next 5-10 years, depending on the weather condition and the 
fraction of fiber time allocated for this survey 0. 



This estimate may still be optimistic. One problem of the Xin- 
glong site is that during the spring to summer time when much 
of the SDSS surveyed region can be observed at meridian transit, 
the weather condition is worse than the autumn (H. Wu, private 
communication) . 



Figure B2. The clustering bias of the main sample as a function 
of redshift. Curves from the top down correspond to the magni- 
tude limit of 17.8, 18.3, 18.8, 19.2, 19.8. 

We can estimate the bias of the main sample using 
the halo model as outlined in §2.3. The HOD parame- 
ter for all galaxy t ypes have been mea sured for the SDSS 
mam sample DR2 l|Zehavi et al.ll2005bl ). which consists of 
204,584 galaxies over 2,947deg^ with apparent luminosity 
10.5 < r < 17.77 and redshifts extended to 2 ~ 0.2. The 
mean number of galaxies above limiting luminosity Lum in 
halos of mass m is parameterized as 

{N[m,L > Lum)) = e(m - mn,in)[l + (m/mi)"], (B6) 

where Q{x) is the Heaviside step function, and the HOD 
parameters 9 = {rrimin , mi , a} are all functions of the lim- 
iting luminosity Lum and re dshift z, their measu red val- 
ues are given m Table 3 of IZehavi et~aLl (|2005bh . How- 
ever, IZehavi et al] |2005b') gives only the measured value at 
Zm ~ 0.1, at higher redshift the data is not available. Since 
the LAMOST galaxy sample would have much higher aver- 
age redshift, evolution of HOD parameters should be consid- 
ered. We have adopted the simple ansatz that the HOD pa- 
rameters are functions of one parameter A only, A is a func- 
tion of both the actual limiting luminosity L and the redshift 
z. At z = Zm, ML, Zm) = L, and 6[X{L,z)] — 6{L,Zm) as 
given m IZehavi et af] (|2005bl ). With this ansatz, in order 
to know 9{L, z) at a higher z, we only need to know the 
corresponding X{L,z). To obtain A, we calculate the galaxy 
density at z using the HOD parameters obtained with a trial 
value of A, compare it with the density integ rated from the 
measured luminosity function (j>{L, z) given in lBlanton et al.l 
(l2003), 

poo J poo 

J dm-£-{N{m,L> Lum,\)) = j (t>{L,z)dL, (B7) 

where {N{m, L > Lum, A)) is given in eq. IB6I The value of 
A is adjusted so that the galaxy densities calculated with 
the HOD matches the observation. Of course, this ansatz 
may not be true. However, since we match the luminosity 
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function at difTerent redshifts, the estimate of bias would not 
be too far off. 

By using this method, we obtain an estimate of the bias 
at different redshifts. For example, for the MAINl sample, 
aX z = 0.1, the bias is 6 = 1.0, and at z = 0.3, b = 1.44. For 
the MAIN2 sample, b = 0.92 at z ^ 0.1, b = 1.14 at z = 0.3, 
and b = 1.66 at z = 0.5. The bias as a function redshift is 
plotted in Fig. IB2I For reference, besides the r = 18.8 and 
19.8 main survey samples, we also plotted several samples 
of different magnitude limits. 

The distributions of galaxies as a function of redshift are 
plotted in Fig.[T] The MAINl sample peaks at z 0.2, but 
extends beyond z = 0.4, where comoving number density 
decrease to n{z = 0.4) ^ 7 x 10"*(/i/Mpc)^. The MAIN2 
sample has a more extended distribution, up to z = 0.6 
with the same cut-off density. The corresponding comoving 
density distributions are plotted in Fig. IB3I 

B2 The LRG survey 

The LRGs are intrinsically brighter and can be selected by 
their colors. At present, a catalogue of LRGs have already 
been compiled from the SPSS photometric da ta, namely the 
MegaZ LRG catalogue (|Collister et. al.ll2006 l). This is a cat- 
alogue of LRG galaxies selected from the SDSS DR4 imaging 
data, which covers 5914 deg'^, with 17.5 < idev < 20.0. The 
total number of galaxies in this catalogue is 1.2 million, and 
photometric redshifts for these galaxies have been obtained, 
with an error of o"z = 0.048. It is expected that the sam- 
ple will be extended to 8000 deg^, and a total number of 
1.5 million of LRGs. The sample selection of LRGs might 
be further tailored for the needs of LAMOST, nevertheless 
it should have a lot in common with the MegaZ sample, 
particularly at 2 > 0.4. We therefore tentatively take the 
properties of this sample as the LAMOST LRG survey sam- 
ple in our estimation. 

We expect that an exposure time of at least 70 minutes 
is required for observing most of these targets. The angular 
density of the sample is about 205 per square degree, which 
is very close to the angular number density of the fibers of 
LAMOST (200). The total observing fiber time required for 
completing 8000 square degree area is about 1.9 million fiber 
hours (80 nights). 

Note that there is not a unique selection criterion of 
LRGs, hence the LRGs in the MegaZ sample may dif- 
fer slightly from those already spectroscopically observed 
in the SDSS-II, or those observed in the 2SLAQ survey 
l|Cannon et al.l [20061 ). The color cut used by MegaZ is sim- 
ilar to but somewhat broader than that used by 2SLAQ. 
It includes bluer objects, which is in turn bluer than the 
Cut II of SD SS LRG to accomm odate the higher density 
of 2dF fibers. IWake et al.l l|2006l l has compared the selec- 
tion criterion difference between SDSS and 2SLAQ for the 
evolution of luminosity function, they started from SDSS 
main sample which consists of all galaxies with r < 17.77, 
then evolved these to 2; = 0.2 (for SDSS) and z — 0.55 
(for 2SLAQ) by interpolating between passive and continous 
star formation model according to the observed g — i color. 
They found that 90% of the SDSS LRG would be selected 
by 2SLAQ, while only 30% of the 2SLAQ LRG would be 
selected by SDSS. The number distribution of the MegaZ 
sample (marked LAMOST LRG) is shown in Fig. [T] most 
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targets in this catalogue have redshifts between 0.4 and 0.7. 
It peaks a roud z = 0.45. Nearly 20% of the objects in the 
catalogue (|Collister et. al.ll2006l ) are located within the re- 
gion with 0.50 < dx < 0.55, d± = {r - i) - {g - r)/8.0. 
If these objects are removed, then the photometric red- 
shift distributions of the MegaZ-LRG catalogue and the 
2S LAQ evaluation set w ould become similar (c.f. Figure. 12 
in lCoUister et. al.l ([2006')). 

To estimate the bias of the Mega-Z sampl e, we use the 
HOD parameters for red galaxies as given in iBrown et al.l 
(j2008h . which is obtained by using the observed luminosity 
function and the luminosity- dependent clustering of 40696 
red galaxies with 0.2 < z < 1.0. The HOD parameters can be 
approximated as a function of redshift and B-band absolute 
magnitude Mb + 1.22, which is an effective proxy for stellar 
mass: 

M„,„(ft_lM0) = 10"-^^ X 10«'4x[-19-(M«-51og;,+1.2.)l 

_|_]^q13.70 ^ -^Ql.l5x[-21-(A/i5-51ogh+1.2z)l _j_ -j^gll.85 j-gg^ 

_(_-]^q14.60 ^ j^Q0.85x[-21-{MB-51ogh+1.2z)l ^gg^ 

with Mq = Mmin, o"iogA/ = 0.3 and a = 1. Since the bands 
are different and the transformation between them are com- 
plicated, again, we calculate the bias by tuning the B band 
magnitude threshold to match the space density, then use 
the formalism presented in §2.3. We find the biases of the 
sample are 1.77 at z=0.475, and 2.29 at z=0.625. 

B3 QSO survey 

Finally, we consider a QSO survey. The Lya forest of the 
QSO can be a powerful probe of the baryon acoustic os- 
cillation (jMcDonald fc Eisenst"einll2007l ). but here the QSO 
itself is used as tracers. In past surveys such as the 2dF and 
SDSS, the low number density of QSOs prevents it from be- 
ing a good tracer target for measuring BAO features. How- 
ever, with LAMOST the comoving number density can be 
greater, making it more useful for s uch purpose. Using the 
2SLAQ sample ([Richards et al.ll2005h . which consists of 5645 
quasars down to g = 21.85 in 105.7 deg^, we obtain a lu- 
minosity function of the quasars. Unlike the Schechter form 
of galaxies luminosity function, the QS O luminosity func- 
tion is modelled as a double power law l|Croom et al ] |2004l : 
[Richards etHlbOOSh 
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redshift 

Figure B3. The comoving number density of galaxy samples. 
^{Mg,z) = ^Q0.4{c« + l)(Mg-MJ) _^ -^Q0.4{/3+l)(Mg-M*) (^^0) 

For estimation we assume here that the evolution with 
redshift is pure luminosity evolution, individual quasars is 
fainter today than at z = 2, with th e dependence of the c har- 
acteristic luminosity described by (iRichards et al.|[2005l ) 

M*{z) = M; - 2.5(fci2 + fc2z^) (Bll) 

where a = -3.28, f3 = -1.78, M* = -22.68, fci = 
1.37, k2 = -0.32, $* = 5.96 x lO'^Mpc-^ In reality, 
the evolution of quasar luminosity function may be due to 
a combination of decreasing space density and decreasing 
luminosity of individual quasars, but this would not signifi- 
cantly affect our result. We calculated the QSO K-correction 
by assuming that the continuum Fi, oc v~'^'^ . 

The quasar number density is plotted in Fig. IB4I Con- 
sider samples with the magnitude limits, g < 20.5 (QSOl), 
g < 21.0 (QS02), and g < 21.65 (QS03), aU between 
z — 0.4 and z = 2.1. The density of these magnitude 
limited samples drops at low redshift, because there is a 
lower limit Mg < —22.5 on the absolute luminosity for 
quasars (this also reflects a decrease in quasar activity at 
low redshifts). The turning redshift zt is the place where 
Mg(zt,mgjim) = Mg^iim, aud Mg,um = -22.5. They are 

Zt,QS01 = 0.75, Zt,QS02 = 0.92, Zt,QS03 = 1-20. 

The angular density of the QSOl is 30deg~^, and the 
total number is 0.24 million for 8000 deg^. The observing 
time is about 15 minutes for signal-to-noise ratio S/N — 4, 
so it requires about 0.06 million fiber hours. For QS02 sam- 
ple, the angular density is 45deg~^, the spectrum can be 
obtained within 30 minutes for the same S/N, and requires 
0.18 million fiber hours. The angular density of QS03 is 
72 deg~^, the required integration time is 1.5 hours, depend- 
ing on the sky background. It needs 0.855 million fiber hours 
to complete the survey. 

The bias of a QSO sample can be estimated using the 
method given in §2.3. Alternatively, it can also be measured 
directly from the pro jected two point correlation function. 
IPorciani et al.l ()2004l ') made an measurement with about 
14,000 quasars at 0.8 < z < 2.1 in the 2dF/6dF QSO 
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Figure B4. Comoving number density of the QSO samples with 
QSOl (g < 20.5), QS02 (g < 21.0) and QS03 (g < 21.65). 

Redshift Survey. They found a slightly greater bias (see Ta- 
ble. |B2} than that of the estimate obtained using the method 
of §2.3. The difference is greater at higher redshifts (1 0% at 
z = 1.5, and 25% at z ^ 1.9) (jMarulh et al.ll2006l ). Our 
estimate may be improved by using more detailed semi- 
analytical models, but for our purpose this simple estimate 
should suffice. 
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